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ABSTRACT 

Long duration gamma-ray bursts are commonly associated with the deaths of mas- 
sive stars. Spectroscopic studies using the afterglow as a light source provide a unique 
opportunity to unveil the medium surrounding it, probing the densest region of their 
galaxies. This material is usually in a low ionisation state and at large distances from 
the burst site, hence representing the normal interstellar medium in the galaxy. Here 
we present the case of GRB 090426 at z = 2.609, whose optical spectrum indicates an 
almost fully ionised medium together with a low column density of neutral hydrogen. 
For the first time, we also observe variations in the Lya absorption line. Photoionisa- 
tion modeling shows that we are probing material from the vicinity of the burst (~ 80 
pc). The host galaxy is a complex of two luminous interacting galaxies, which might 
suggest that this burst could have occurred in an isolated star-forming region outside 
its host galaxy created in the interaction of the two galaxies. 

Key words: gamma-rays: bursts - X-rays: general - X-ray: ISM 



1 INTRODUCTION 

Optical afterglow studies of gamma-ray bursts (GRBs) pro- 
vide a powerful tool for unveiling the interstellar medium 
(ISM) properties of their host galaxies. Long GRBs are 
connected to the death of a massive star, and their host 
galaxies are subsequently expected to be sites of heavy star- 
formation extending back to the early Universe. Metallicities 
determined from absorption lines of the ISM show higher 
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values than those of, for example, QSO absorbers (Fynbo et 
al. 2006a; Savaglio et al. 2006; Savaglio 2009), though this 
might be only a sightline effect (Fynbo et al. 2008). Kine- 
matics of absorption lines give some indications for galactic 
outflows (e.g. Thone et al. 2007), which are an expected 
phenomenon in star-forming galaxies. 

Usually, absorption lines in GRB afterglow spectra are 
constant in time. In a few cases metal absorption lines to- 
gether with their corresponding fine structure transitions, 
caused by UV pumping from the intense afterglow radia- 
tion, vary in intensity following the change in the ionizing 
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flux of the burst. Modeling of the line variability allows us to 
better constrain the distance from the GRB to the absorb- 
ing material, which typically ranges from > 150 pc to a few 
kpc (e.g. Vreeswijk et al. 2007, D'Elia et al. 2009). The pres- 
ence of neutral material in the spectra implies distances on 
the order of a few hundred pc to several kpc (e.g. Prochaska 
et al. 2007, Thone et al. 2007). For strong NV absorption 
lines, it has been suggested that those lines must arise from 
material close to the GRB site, at <100pc (Prochaska et 
al. 2007). GRB 080310 showed variation in all resonant Fe 
II transitions which can only be explained with photoion- 
ization by the burst at a distance of ~ 100 pc (A. De Cia 
et al. in prep.). Excluding those few special cases, optical 
GRB afterglow spectroscopy seems to predominantly probe 
the ISM of the host galaxy and not the material around the 
GRB progenitor. 

Studying the circumburst environment, however, would 
allow us to verify the predictions of current progenitor mod- 
els. Prompt emission X-ray data showed some evidence for 
variations in the absorbing column density (Amati et al. 
2000). In the optical, absorption features from circumburst 
material are difficult to detect due to the flash ionization of 
the surrounding material by the GRB out to several tens 
of parsecs (Robinson et al. 2010), even though the GRB is 
predicted to reside in a dense environment. In addition, any 
remaining absorbing material would be hidden under the 
stronger absorption lines from the line-of-sight ISM in the 
host galaxy. 

There are several potential signatures which could 
demonstrate that we are observing material close to the 
GRB instead of the general ISM of the galaxy: time variabil- 
ity of the absorption line strength (Perna & Loeb 1998), a 
high ionization state of the medium (Prochaska et al. 2008), 
or absorption lines at large velocities compared to the red- 
shift of the GRB. The last has been observed in a few cases 
(Fox et al. 2008, M0ller et al. 2002), but it cannot be ex- 
cluded that those systems are intervening absorbers and 
hence not connected to the GRB. It was also suggested that 
a large difference between the optical and X-ray column den- 
sities might be a hint of observing highly ionized material, 
eventually close to the GRB (Watson et al. 2007, Campana 
et al. 2010). An ideal laboratory would be a GRB progeni- 
tor that exploded in a relatively isolated environment with 
little ISM along the line of sight. So far, no conclusive evi- 
dence for such conditions has been found, probably because 
such physical conditions are rare and only a few afterglows 
have time resolved afterglow spectroscopy. In this paper we 
present the first example for such a scenario, GRB 090426. 

GRB 090426 was discovered by the BAT (Burst 
Alert Telescope) 7-ray telescope onboard the Swift satellite 
(Gehrels et al. 2004) on AprU 26, 2009, 12:48 UT (Cummings 
et al. 2009) and had a duration of Tgo = 1.2 ± 0.3 s (Sato 
et al. 2009). Soon after its discovery an optical counterpart 
was found by Swift (Cummings et al. 2009) and the redshift 
was determined from an optical spectrum to be 2: = 2.609 
(Levesque et al. 2010). Both its restframe duration (0.3s) 
as well as its observed duration put GRB 090426 into the 
short burst category according to the commonly used clas- 
sification based on the duration of the prompt emission. 

In §2 we describe the optical spectra and afterglow 
observations. In §3 we discuss the characteristics of GRB 
090426 itself and some issues regarding its classification. §4 



describes the properties of the ISM in the line of sight as de- 
rived from optical and X-ray observations. In §5 we deal in 
detail with the Lya line variation and model its behaviour 
with a photoionisation code. §6 describes the properties of 
the host galaxy of GRB 090426. Discussion and conclusions 
are reported in §7. 



2 OBSERVATIONS 

We obtained both photometric and spectroscopic observa- 
tions of the afterglow and the host galaxy. The imaging data 
are listed in Tab. IT] and the spectroscopic observations are 
detailed in Tab. |2| We observed the optical afterglow with 
MOSCA and AlFOSC at the NOT as weU as F0RS2 at the 
VLT from 2009 Apr. 26 to June 19 in several bands. In ad- 
dition, we include acquisition images from the spectroscopic 
observations. Photometry was performed using SDSS stars 
(Abazajian et al. 2009) as reference and the transformations 
described by Jester et al. (2005) to convert the SDSS to V 
and Rc magnitudes. In order to account for the contribution 
of the host galaxy system to the light curve independent of 
the seeing, we used aperture photometry centered on the 
afterglow with radius of 2" . 

We obtained a spectrum of the afterglow with F0RS2 
at the VLT starting on 2009 Apr 27.048 UT (12.3 hr after 
the burst) using grisms 600V and 600RI together with a I'.'O 
slit and an exposure time of 30min per grism. The spec- 
trum covers a wavelength range from 3500 to 8600 A with a 
resolution of ~ 11 A. Furthermore, we reanalyzed the spec- 
trum from Levesque et al. (2010) taken with LRIS at Keck 
1.1 hr after the burst. This spectrum has a similar resolu- 
tion and wavelength coverage compared to the FORS2 spec- 
trum. The spectra were reduced using standard tasks 
in IRAF, with the different exposures and grisms 
from each instrument combined, flux calibrated, and 
then normalized. Flux calibration prior to normal- 
izing the spectrum makes the normahzation more 
reliable, especially bluewards of Ly a and in the 
Ly a forest, since the afterglow continuum follows 
a smooth power law. 

After the fading of the afterglow an object was found 
I'.'l SW of the GRB position and initially suggested as the 
candidate host galaxy (Gl). We took spectra of Gl with 
FORS2 at the VLT on June 18 2009 using grism 300V (3300 
- 8700 A) and exposure times of 2 x 1800 s. The acquisition 
image for the spectroscopic observations of Gl revealed an- 
other object underlying the GRB position (therefore sug- 
gested as the real host galaxy, H), for which we obtained 
another spectrum on June 19 2009 with the same settings 
and exposure times. We also obtained 2 x 30min spectra 
with X-shooter at the VLT (D'Odorico et al. 2006) on June 
2 2009, covering a wavelength range from 9800 to 24800 A. 
The slit was placed across the GRB position (and the host) 
and Gl. 



3 GRB 090426: SHORT OR LONG? 

The original division of GRBs into two populations dates 
back to Kouveliotou et al. (1993) who found a bimodal 
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Table 1. Log of the imaging observations. 



Instrument 


Start time (UT) 


Time after GRB (d) 


Exp. time (s) 


Band 


Magnitude 


MOSCA/NOT 


2009 


Apr 26.879 


0.345 


3x300 


R 


21.27±0.16 


MOSCA/NOT 


2009 


Apr 26.939 


0.405 


1x300 


R 


21.33±0.11 


MOSCA/NOT 


2009 


Apr 26.946 


0.412 


2x300 


i 


21.18±0.02 


MOSCA/NOT 


2009 


Apr 26.958 


0.424 


5x200 


z 


21.12±0.07 


MOSCA/NOT 


2009 


Apr 27.015 


0.481 


2x300 


g 


21.86±0.03 


FORS2/VLT 


2009 


Apr 27.048 


0.514 


1x120 


R 


21.52±0.04 


MOSCA/NOT 


2009 


Apr 27.113 


0.579 


2x300 


g 


22.12±0.05 


MOSCA/NO i 


2009 


Apr 27.122 


0.588 


2x300 


V 


22. 05±0. 06 


MOSCA/NOT 


2009 


Apr 27.128 


0.594 


1x300 


R 


21.77±0.12 


FORS2/VLT 


2009 


Apr 27.144 


0.610 


1x120 


R 


21.92±0.04 


FORS2/VLT 


2009 


Apr 27.150 


0.616 


1x120 


R 


21.95±0.04 


FORS2/VLT 


2009 


Apr 27.152 


0.618 


1x60 


R 


21.87±0.04 


MOSCA/NOT 


2009 


Apr 27.883 


1.349 


3x300 


R 


22.85±0.19 


MOSCA/NOT 


2009 


Apr 28.159 


1.625 


4x300 


R 


23.00±0.13 


MOSCA/NOT 


2009 


Apr 29.067 


2.533 


6x600 


R 


23.30±0.11 


AlFOSC/NOT 


2009 


May 03.918 


7.384 


10x600 


R 


23.73±0.19 


FORS2/VLT 


2009 


Jun 18.981 


53.447 


2x300 


V 


23.76±0.17 


AlFOSC/NOT 


2009 


Jun 19.915 


54.381 


7x600 


R 


23.71±0.17 


FORS2/VLT 


2009 


Jun 19.989 


54.455 


300+450 


V 


23.88±0.17 


FORS2/VLT 


2009 


Jun 19.976 


54.442 


1x300 


R 


23.71±0.31 



Table 2. Log of the spectroscopic observations. 



Instrument 


Start time (UT) 


Exp. time (s) 


Grism 


Resolution 


Keck/LRIS 


2009 Apr 26.506 


2x300 


300/5000 


650 


VLT/F0RS2 


2009 Apr 27.056 


2x1800 


600B 


740 


VLT/FORS2 


2009 Apr 27.100 


2x1800 


600RI 


740 


VLT/X-shooter 


2009 May 31.975 


3x1200 


VIS+NIR 


~6000 


VLT/FORS2 


2009 Jun 18.997 


2x1800 


300V 


440 


VLT/FORS2 


2009 Jun 20.011 


2x1800 


300V 


440 



distribution in the high-energy duration and hardness ra- 
tio (HR) of BATSE bursts. This was thought to imply two 
different progenitor populations, in particular since collap- 
sars cannot produce bursts with a duration shorter than a 
few seconds. A third, intermediate, population was also pro- 
posed on statistical grounds (Horvath et al. 2010), but there 
is no conclusive evidence from properties of the bursts, af- 
terglows, or host galaxies that this third class is intrinsically 
different from long bursts (de Ugarte Postigo et al. 2010). In 
recent years, the classification of GRBs based only on their 
high-energy properties has been widely discussed after the 
discovery of outliers such as long GRBs without detected 
supernovae (SNe) and a number of high-redshift bursts that 
are short in their rest-frame duration but likely coming from 
the collapse of a massive star. Levesque et al. (2010) statis- 
tically placed GRB 090426 in the short burst category; how- 
ever, due to other properties they did not exclude a coUapsar 
progenitor model. 

In the following, we present some properties of the burst 
and afterglow that suggest it was likely produced by a coUap- 
sar event. In order to compare GRB 090426 with the BATSE 
bursts we transform the count ratio of the Swift-YiNT data 
to those of the BATSE channels and obtain a ratio of (50- 
100 keV) / (15-25 keV) = 0.60, which classifies it as a "soft" 
burst. This places GRB 090426 in a very exceptional region 
of the HR vs. Tgo diagram, outside of the two main groups of 
long-soft and short-hard bursts with a spectrum that is even 



softer than most long-duration Swift bursts (see Fig.[T|. The 
luminosity of its prompt emission with E\ao = (5± 1) x 10^^ 
erg (1-10 000 keV, isotropic energy release, Antonelli et al. 
2009) would be considered rather high for a merger event. 
However, it could be possible to reach such high 
luminosities through merger events. Recent simula- 
tions of neutron star-black hole mergers (Pannarale, 
Tonita & Rezzolla 2010) show that the torus created 
in the merging process can approach or even exceed 
one solar mass, which would yield enough energy to 
power this burst. A high energy release could also 
be facilitated by neutrino-powered outflows or mag- 
netic fields (e.g. Aloy, Janka & Miiller, 2005). 

Another distinction between long and short burst is 
the determination of a spectral lag between the high energy 
emission in different bands, which is usually present for long 
bursts but absent for short bursts (Norris et al. 2006). For 
GRB 090426, the lag is consistent with zero but has a large 
uncertainty (Ukwatta et al. 2009). 

The lightcurve of the optical afterglow was previously 
discussed in AntoneUi et al. (2009) and Xin et al. (2011). We 
collect the lightcurve data available from the literature from 
90 s to 35 d after the burst and add the observations listed 
in Tab. [l] We normalize all optical data to the ii-band to 
construct a single lightcurve. The optical light cuve of GRB 
090426 is similar to what is usually found for other long (or 
short) GRB afterglows. It reaches its maximum before 100 s. 
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Figure 1. Tqq vs. hardness-ratio according to the BATSE chan- 
nels (50-100 keV vs. 15-25 keV) applied to Swift bursts with 
known redshift (Horvath et al. 2010). Bursts marked in black 
are those with a column density of neutral hydrogen determined 
from optical spectra with log Njji/cm"^ < 20. 

when the observations pubhshed by Xin et al. (2011) begin, 
which likely marks the onset of the afterglow. Between 300 
and 3,000 s the light curve decays with a temporal slope of 
1.12 ± 0.02. At 3,000s we see a flattening in the light curve 
after which the decay seems to continue with approximately 
the same slope as before. This bump could be due to an 
energy injection (Bjornsson et al. 2004), as it is achromatic 
and happens simultaneously in optical and in X-rays (see 
Fig.[2]|. Similar behaviour has been seen in a number of other 
bursts e.g. GRB 021004 (de Ugarte Postigo et al. 2004) and 
GRB 060526 (Thone et al. 2010). There are also indications 
of similar bumps at 10,000 and 30,000 s. After the first day, 
the lightcurve flattens when the liglit of the host galaxy at 
7? ~ 23.8 mag begins to dominate. In Fig. [2] we show the 
optical together with the X-ray lightcurve in flux units. The 
optical data were not corrected for the contribution from the 
host galaxy. 

Kann et al. (2010) found a clear bimodality in the dis- 
tribution of optical afterglow luminosities, with an average 
absolute magnitude of Mb= -23.45 mag for long/type II 
GRBs and Mb= -18.2 mag for short/type I GRBs at 1 day 
after the burst (restframe). This might be due to dif- 
ferences in the density of the surrounding material, 
which should be higher for long bursts that are typ- 
ically found in dense star-forming regions. With an 
absolute magnitude of Mb= —23.79 mag, the lumi- 
nosity of the afterglow of GRB 090426 is around 
the average for long GRB afterglows, and about 5 
magnitudes brighter than a typical short GRB. It 
has also been suggested that some compact object 



Figure 2. Optical lightcurve collected from data in the litera- 
ture (Antonelli et al. 2009; Xin et al. 2011; open symbols) and 
observations from the NOT and VLT (published in this article; 
filled symbols), all normalized to the R-band. The dashed line 
indicates the flux contribution from the underlying host galaxy. 
We plot the X-ray lightcurve for comparison with the flux density 
provided by the Swift archive (Evans et al. 2010). 



mergers have a very short coalescence time while 
they are still within their star forming region, and 
hence should have afterglow luminosities compara- 
ble to long bursts (Belczynski et al. 2006). The rela- 
tively high redshift of this burst compared to other 
short bursts would support this possibility. 

As previously suggested by other authors, we consider 
this burst to be peculiar. However, the properties detailed 
above place it closer to the population of GRBs originating 
from massive progenitors. 



4 PROPERTIES OF THE ISM IN THE 
LINE-OF-SIGHT 

4.1 X— ray spectra analysis 

We analysed seven observations with the XRT instrument 
onboard the Swift satellite, covering a time period from 26 
April to 5 May 2009. The source was only detected during 
the flrst two observations on 26 April 2009. The XRT data 
were reprocessed with standard procedures (xrtpipeline ver. 
0.12.4 within the Heasoft package ver. 6.8), resulting in 4.7 
and 2.0 ks exposures in photon counting mode, respectively. 
We decided to group the events to one photon per channel 
and to adopt Cash-statistics (Cash 1979) within the fltting 
program XSPEC. At the beginning of the first observation 
the source is slightly piled-up (two photons fall on the same 
pixel during the same integration time and are read as a 
single photon with higher energy). We take a conservative 
pile-up threshold of 0.7 ct s~^ and split the first observation 
into a piled-up interval (1.8 ks) and a non-piled-up interval 
(2.9 ks) which was summed to the second observation for 
a total time of 4.9 ks. Photons from the piled-up interval 
were extracted from an annular region with 2 and 25 pixel 
inner and outer radii centered on source. Photons from the 
non-piled-up observations were extracted from a 20 pixel 
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radius circular region. The background was extracted from 
the same annular region with an inner and outer radius of 
60 and 80 pixels, respectively, both of which are free of field 
sources. 

Ancillary response files were generated with the 
xrtmkarf task, accounting for CCD defects, vignetting, and 
point-spread function corrections using the appropriate ex- 
posure map. We end up with two spectral files containing 
383 and 328 counts in the 0.3-10 keV energy band. The back- 
ground is very low, comprising of less than 1% and 7% for 
the two observations, respectively. We group the events and 
adopt Cash-statistics (Cash 1979) within the fitting program 
XSPEC. During the spectral fit we used the newly released 
broadened response file swxpc0tol2s6_20070901v011 .rmf . 
The spectral model consists of a power law model and 
two absorption components, one from our Galaxy fixed to 
Nh = 1.5 X 10^" cm"^ (Kalberla et al. 2005) and one at the 
redshift of the host galaxy {z = 2.609). 

We find consistent spectral results between the two 
spectra, so we fit the column densities and the power law 
photon indices together and leave only the normalizations 
free to vary. Assuming solar metallicity and solar compo- 
sition, the best fit column density is Nh ~ 5.0t-xl x 10^^ 
cm~^ and the power law photon index is F = 2.1 ± 0.1 (er- 
rors are at 90% confidence level with Ax^ = 2.71)s. The 
column density is around the average found for other long 
Sunft GRBs (Campana et al. 2010). 



4.2 The optical spectrum 

The spectrum of the optical afterglow turns out to be very 
different from most other GRB spectra. In our FORS spec- 
trum we only detect very strong high ionization lines of N V 
AA 1238, 1242, Si IV AA 1393, 1402, and CIV AA 1548, 1550. 
In addition, the Lya absorption is relatively weak; however, 
this is not unprecedented in GRB afterglow spectra. In con- 
trast to the usually strong low ionization lines detected in 
other GRB spectra, we only detect a weak CII A 1334 line 
while normally prominent lines of Si II AA 1260, 1526, Fell 
and AlII are undetected. This implies that the observed 
medium is in a high ionization state, which is very unusual. 

We then compared our FORS spectrum from 12 h after 
the burst with the the spectra by Levesque et al. (2010) from 
1.1 h after the burst. We detect the same high ionization 
absorption lines, but the Lya line is significantly stronger 
than in our later spectra. We measure a 3.8 a decrease in 
the equivalent width (EW) between the two spectra at 1.1 
and 12 h after the GRB (see Fig. [sj. The other absorption 
lines do not show any significant variability. Variability of 
Lyman a has never been before observed in a GRB afterglow 
spectrum, nor in any other extragalactic object. 

The column densities were obtained from the curve-of- 
growth (CoGQ adopting a 6-parameter of 50 kms"'^. The 
errors on the Lya column densities only refiect the measured 
errors on the EWs, not the uncertainty in the 6-parameter, 



^ The CoG describes the evolution of the column density with 
EW, with a linear relation for unsaturated lines and a dependence 
on the broadening of the line by turbulence or temperature in the 
following, flat part of the curve. 
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Figure 3. Top panel: Lya absorption lino at two different epochs; 
the dotted lines are the corresponding error spectra. Bottom 
panel: Difi^erence between the two epochs (solid line) and 1 cr 
errors (dotted lines). The spectra have been smoothed with a 
running window before subtraction, corresponding to the instru- 
mental width. 



which would change the absolute values for the column den- 
sities but not the significance of the temporal variation. For 
the other lines, we assume very conservative errors on the 
column densities to account for the uncertainty on the b- 
parameter of around 15 kms~^. Si IV has two components 
resolved and we therefore determine the total column den- 
sity by summing the contributions from the individual com- 
ponents (which are each closer to the linear part of the 
CoG than the combined line). For CIV, the column density 
was calculated assuming an EW ratio of 1:2 for the blended 
1549/1550 doublet and taking the column density from the 
CoG for the weaker transition (1550 A). The column density 
listed is therefore not for the doublet together but closer to 
the "true" column density (had it been possible to deblend 
the two lines). However, since the line is clearly saturated, 
the true ratio deviates from the theoretical 1:2 ratio and the 
error in the column density must therefore be considered 
much larger. 

In Table 3 we list the measured equivalent widths 
(EWs) and derived column densities for the VLT/FORS 
spectrum and the comparison with the spectrum presented 
in Levesque et al. (2010). We also list the deviation of the 
individual lines between the two epochs. Both spectra have 
been analysed consistently by our group, and the EW values 
might therefore slightly deviate from those listed in Levesque 
et al. (2010). A plot of the F0RS2 spectrum is shown in Fig. 

H 

The variability of Lya strongly suggests that the mate- 
rial lies close to the GRB progenitor. Such an effect has been 
predicted from photoionization modeling (Perna & Loeb 
1998), but has never before been observed. The high column 
densities of the high ionization lines and the absence of their 
low ionization counterparts support the scenario that we are 
observing a large quantity of matter close to and ionised by 
the GRB. This is further supported by the large difference 
between the optical and X-ray column density, which can 
be explained by photoionization of the medium by the GRB; 
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Observed wavelength (A) 

Figure 4. Normalized spectrum from FORS2/VLT grism 600B and 600RI with absorption lines indicated. Undetected lines that are 
commonly detected in other GRB spectra are written in brackets. Atm. A and B mark the major atmospheric absorption bands. 



Table 3. Absorption line properties and upper limits on some transitions that arc usually strong in other 
GRB spectra but undetected here. EWs are given in restframe. 



Ion 


■^rcst 

(A) 


'^obs 

(A) 


EWvLT (12 h) 
(A) 


log N 
(cm"^) 


(A) 


(l.lh) 


EW diff. 


HI 


1215 


4385 


1.94±0.23 


18.7lO-i 


3.09±0.2 




-3.8(7 










(log Nhi = 


19.2 ±0.1) 




NV 


1238 


4471 


0.99±0.20 


15.8±0.8 


0.70±0.10 






NV 


1242 


4485 


0.29±0.07 


15.8±0.8 


0.30±0.10 






Sill/SII 


1260 


5457 


< 0.35 


< 13.6 


< 0.35 






Sill 


1304 


4797 


< 0.25 


< 13.6 








CII 


1334 


4814 


0.86±0.19 


15.8±1.0 


0.53±0.14 




+ 1.3(T 


Si IV 


1393 


5030 


2.09±0.20 


15.4±1.0 


1.78±0.15 






Si IV 


1403 


5062 


1.38±0.20 


15.4±1.0 


1.53±0.14 




-1.5o- 


Sill 


1526 


5504 


< 0.35 


< 13.6 








CIV 


1549 


5591 


4.32±0.28 


17.6±1.0 


3.63±0.25 




+ 1.8(7 


Fell 


1608 


5804 


< 0.34 


< 14.4 








AlII 


1670 


6029 


< 0.37 


< 12.9 









X-ray absorption is insensitive to this phenomenon since it 
is caused by excitation of the electrons in the inner shells of 
the atoms (Campana et al. 2010). 



4-2.1 Comparison to the "average" GRB spectrum 

We compare the lines found in the spectrum of GRB 090426 
to a sample of 77 low-resolution afterglow spectra of long 
GRBs observed by Swift (Fynbo et al. 2009). The majority 
of those spectra in the right wavelength range to detect Lya 
absorption (33) showed a strong damped Lya (DLA) system 
with a median value of log Nni/cm^^ = 21.6. Only 6 other 
spectra had column densities below log Nni/cm^^ — 20.0, 
the lower limit for the definition of a DLA: GRB 050908 (log 
Nm/cm-^ = 17.8), GRB 060124 (log Nm/cm-^ = 18.5), 
GRB 060607A (log Nni/cm-^ ~ 17.2), GRB 070411 (log 
NHi/cm-2 = 19.3), GRB 071020 (log Nm/cm-^ < 20.0), 
and GRB 080310 (log Nni/cm"^ = 18.8). All of them se- 
curely belong to the long burst group with Tgos between 25 
and 100 s. Five of these bursts are also highlighted in Fig.[l] 
Christensen et al. (2010) produced a composite using 60 
of the 77 spectra from Fynbo et al. (2009) mentioned above, 
and determined EWs and column densities for all detected 



lines in the composite. The following lines, which are usually 
strong in the environments of long-duration GRBs, are not 
detected in our spectrum: Sill A 1260,1304 and 1526 A, C 
II A1534, and Al III A 1670,1855 and 1863. The EWs of C 
IV and Si IV in the spectrum of GRB 090426 are almost 
twice as high as in the composite GRB spectrum. Consider- 
ing that the HI column density of GRB 090426 is very low 
compared to the approximate value for the composite (log 
Ai'Hi/cm~^= 21.7), the relative abundance of these species 
are rather extreme for this burst. 

If we take only a composite of sub-DLA bursts (3 spec- 
tra), we see the following: a fit of the Lya absorption line 
gives log A''Hi/cm~'^=19.3±0.1 (dominated by the spectrum 
of GRB 070411). All the usual strong low ionization lines are 
also present in this composite such as Si II A 1260, 1304 and 
1526 A, Si II* A 1264 A, C II A 1334 A, and Al III A 1670 
A. Si IV A 1402 is not very strong, but falls in a noisy area 
of the composite, with EWrost= 0.25±0.20 A. In contrast, 
CIV 1549 is also quite strong here, with EWrost of 3.2±0.1 
A. The values in Fig. [H] are not derived from the composite 
spectrum since this is dominated by the GRB 070411 spec- 
trum; instead instead we took the average of the EWs of all 
sub-DLAs. 
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Comparing the ionisation of C and Si to those deter- 
mined from the sample of Fynbo et al. 2009 (see Fig. [5|, 
GRB 090426 occupies a rather different region than most 
burst in the sample. However, it seems that other GRBs 
with low Nhi also tend to show a rather high ionisation state 
of the observed medium. GRB 060607A is the most extreme 
case, although it is not included in the plot since only CIV 
has been listed in Fynbo et al. (2009). From Fig. 3 in Fox et 
al. (2008), who analysed the UVES spectrum of this burst, 
we see that this burst, together with the lowest reported 
log Nni/cm"^ of 16.8 (Fox et al. 2008, whereas Fynbo et 
al. 2009 report a value of ~ 17.2) for a GRB, has no CII 
or Sill while the high ionisation species are rather strong. 
GRB 050809 has relatively strong low ionisation lines, but 
the value in Fig. [5] is influenced by the nondetection of Si IV 
A 1403 while the other line of the doublet at A 1393 is de- 
tected. This seems to confirm that there is a correlation 
between the hydrogen column density and the ionisation of 
the medium, not surprising since a large hydrogen column 
density would screen the UV radiation from the GRB. How- 
ever, it might also imply that some bursts occur in a rather 
isolated region of their host galaxy, which would allow us to 
truly probe ionised gas from the surroundings of the burst 
itself. 



5 PHOTOIONIZATION MODELING 

Photoionization models allow to us determine the distance 
of the observed material from the GRB by considering the 
ionization state, column densities, and their variations with 
time (e.g. Lazzati et al. 2001). Here, we use the photoion- 
isation code by Perna & Lazzati (2002, see also Lazzati & 
Perna 2003) specifically written for application to GRB ob- 
servations. The code assumes that the GRB is initially sur- 
rounded by a cold (To = 100 K) absorbing material with 
solar metallicity, either in the form of a spherical shell with 
a thickness of 0.1 the shell radius or a spherical cloud. 
The GRB spectrum and its temporal evolution are the input 
in the code, which computes the evolution of the absorbing 
material's ionization state as a function of both time and the 
distance of the absorber from the burst. At each time, the 
resulting wavelength-dependent opacity integrated over the 
line-of-sight is outputted and compared to the observations. 

The total absorption at early times is a measure of the 
total absorbing column, while the speed of the evolution is a 
measure of the distance of the absorber. This analysis is sim- 
ilar to what has been done for the combined soft-X-ray pho- 
toionization opacity in GRB 050904 (Campana et al. 2007). 
GRB 090426 is better suited for this kind of analysis since 
we observe variation in the HI absorption. Since hydrogen 
is the most abundant element in the ISM, constraining the 
geometry of the absorption region with the evolution of the 
Lya absorption gives us a direct result that does not depend 
on metallicity considerations. The four panels of Figure [6] 
show the measured column densities of the H, C, Si, and N 
ions with the best absorption model overlaid. We obtain the 
best fit model for a uniform cloud surrounding the GRB. 
The cloud properties are fit using only the HI absorption, 
with the other ions and the soft-X-ray continuum absorption 
only used as consistency checks. 

The absorber that best reproduces our Lya measure- 
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Figure 5. Top: Comparison of EWs of the average optically 
bright long-duration Swift burst, a subsample of bursts with a 
low neutral hydrogen column density, and the values for GRB 
090426. Also plotted are the ionization ratios of Si and C. This 
comparison clearly shows the exceptional line strengths of the 
spectrum of GRB 090426, showing a highly ionised medium. Bot- 
tom: Ionization of Si vs. C for GRB 090426 and the GRB sample 
presented in Fynbo et al. (2009). Listed are only those bursts 
with a detection in at least one of the ionization states of Si and 
C. Bursts with log Nni/cm-^ < 19 are labeled. GRB 090426 has 
a high ionization for both elements, similar to other bursts with 
low log Nhi- 



ments is a cloud of radius i? = 79 pc (which corresponds to 
the size of a large molecular cloud) and an initial HI column 
density A^hi = 2 x lO^'' cm^"^. Such values correspond to a 
total mass of the molecular cloud of ~ 40, 000 Mq and an 
average density n ~ 1 cm~^. The errors on the values from 
the modeling are on the order of 20%. The model does not 
rule out a solution in the form of a shell; however, 
since the thickness of the shell is O.IR, the density 
would be about 10 times larger as a result. A wind 
blown medium would be fully ionised (see Robinson 
et al. 2010) at very early times and is therefore ruled 
out here. 

Figure [7| shows the 1-, 2-, and 3- a contour levels of the 
Lya fitting in the R — Nm plane (thin solid lines). Overlaid 
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Figure 6. Temporal evolution of the absorbing columns of, from 
top to bottom, H, C, Si, and N ions. The best fit evolution of the 
absorbing columns as derived from our photoionisation model is 
overlaid. 



is the area that is consistent with the combined constraint 
of Si IV column evolution and soft X-ray continuum absorp- 
tion. The Si IV constraint allows for regions of high column 
density (right of the thick solid line) while the soft X-ray ab- 
sorption allows only for low column density regions (left of 
the dashed line). The shaded area is in full agreement with 
the Lya result. A much looser constraint comes from the NV 
fitting, while the CIV column density cannot be fit by our 
model. The C IV excess might be an indication of a carbon- 
enriched medium, possibly from earlier SN explosions, but 
we caution the reader due to the large errors affecting the 
CIV line (see Sec.|4!2|. 

The agreement of the photoionization model predictions 
with the evolution of the column densities ultimately con- 
firms that all the material we observe is close to the GRB 
site. Such material has never been probed before in GRB 
afterglow spectra. 




Log(lnitial column density) (cm ^) 

Figure 7. Constraints in the R — Nm plane from the photoion- 
ization modeling. The thin solid lines show the 1-, 2-, and 3-cr 
contour levels from the Lyo fitting only. The gray area shows 
the region allowed by the combined constraint on the Si IV evo- 
lution (right of the solid thick line, allowing only high column 
density regions) and the continuum soft X-ray absorption (left of 
the dashed line, allowing only low column density regions). 




Figure 8. Color image of the host galaxy field and magnified 
complex of the host galaxy (H) and the companion (Gl). The 
position of the afterglow is marked with a red cross in the inset. 



6 THE HOST GALAXY COMPLEX 

In addition to the afterglow spectroscopy, we obtained ob- 
servations of the host galaxy and its neighbour. Using the 
acquisition images of the spectroscopic observations, we de- 
termined the magnitudes of the two galaxies Gl and H 
(the host) to be Vbi = 24.13 ± 0.05 and Vk = 24.93 ± 
0.08, respectively, corresponding to absolute magnitudes of 
Mab(1500 A)= -21.05 and -20.43. At the host redshift, this 
implies relative luminosities of 1.5 and 0.8 L, (Reddy & 
Steidel 2009). Fig. [s] shows a color image of the host galaxy 
complex (obtained from our acquisition images) and indi- 
cates the position of the afterglow. 

In the F0RS2 observations, both galaxies show clear 
LyQ emission at redshifts of z = 2.605 (H) and 2.611 (Gl), 
indicating that the host has about the same redshift as the 
GRB while the companion is at a slightly higher redshift. 
At that redshift, the physical distance between their cores is 
only 9kpc in projection, which implies that they are likely 
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interacting and possibly in tiie process of merging. The Lya 
emission of the host galaxy gives a star formation rate of 3.1 
Mq yT~^, on the lower end of what is typically observed in 
(long duration) GRB hosts (Savaglio et al. 2009). 

In the 2-dimensional X-shooter spectrum we detect sev- 
eral faint emission lines in the near-IR, including [OIII] A 
5007, 4959 and H/3. Both [O III] lines show two components 
separated by about 400 km s~^ and offset in the spatial di- 
rection, which we attribute to the host (H) and Gl. The 
second component of H/3 is coincident with an atmospheric 
emission line. 

Comparing the afterglow spectra and the emission lines 
from the hot gas of the two galaxies, the kinematics turn 
out to be rather complex (see Fig. [9|. The Si IV transitions 
at A 1393 and 1403 A each show two equally strong velocity 
components with a separation of about 300 kms~^. For N V 
and C IV, the two components cannot be separated since 
the doublets are closer together (C IV is also affected by a 
sky line) . The red component of the Si IV lines is at about 
the same velocity as the center of the Lya absorption. The 
almost equally strong blue Si IV component is coincident 
with the blue wing of Lya, which could be fitted with a 
small additional component at that velocity but with only 
1/100 to 1/1000 the column density of the main compo- 
nent. The Lya emission lines of the host and Gl are slightly 
redshifted compared to their nebular emission lines, which 
could imply an outflow of hydrogen gas as commonly ob- 
served in, for example, Lyman break galaxies (Steidel et al. 
2010). Furthermore, emission and absorption line velocities 
from the host do not coincide, which suggests a complex ve- 
locity held due to the physical motion of the material seen 
in absorption. This is to be expected in the merging process 
of two galaxies. It also supports the suggestion that the ma- 
terial we see in absorption is actually close to the GRB and 
not only coming from the ISM of the host galaxy. 



1^ 






* Host 










Companion 


f 




• 



5,0 
4,0 
3,0 
2,0 
1,0 
100 




-3000 -2000 -1000 1000 2000 3000 

Velocity (km/s) 

Figure 9. Kinematics of the different absorption and emission 
lines from the GRB afterglow spectrum and the host galaxy com- 
plex. Panel (a): Lya in absorption from the FORS2 spectra (12 
hr after the burst) together with the fit and the la errors of the 
Lya profile. Panel (b): Si IV in absorption with two velocity com- 
ponents indicated by the dashed-dotted lines. Panels (c) &; (d): 
Lya in emission from the host and its companion (FORS2). Panel 

(e) : 2D spectrum from X-shooter with the [O III] emission lines 
detected, tracing the ISM in the host and its companion. Panel 

(f) ; an image of the host and its companion with the slit position 
indicated. 



7 DISCUSSION AND CONCLUSIONS 

The spectrum of GRB 090426 is the first compelling case 
of observed photoionised material from the star-forming re- 
gion of the progenitor itself. Its low hydrogen column den- 
sity is varying, the absorption lines show a highly ionised 
medium of high column density, and the large discrepancy 
between X-ray and optical absorption gives another indica- 
tion for photoionised material being present. GRB spectra 
have shown a large range of hydrogen column densities, usu- 
ally explained by different sightlines of the GRB through its 
host galaxy (Fynbo et al. 2009; Jakobsson et al. 2006), with 
some column densities even lower than the one measured for 
GRB 090426. However, if GRB 090426 exploded in, for ex- 
ample, the halo of the galaxy or a very low density ISM like 
a globular cluster, we would also expect low column den- 
sities of the metal absorption lines. A location in the halo 
had been suggested for GRB 070125 (Cenko et al. 2008) 
and GRB 071003 (Perley et al. 2008), both of which show 
a bright but featureless continuum with extremely weak ab- 
sorption lines. GRBs with log TVni/cm"^ < 20.0 tend to 
show a higher fraction of ionised versus neutral material, 
possibly explained by a lack of UV shielding by the lower 
density of hydrogen; however, some of them still contain a 
relatively large column density of neutral material. In ad- 



dition, the afterglow showed that the GRB must have ex- 
ploded in a non-low-density medium (Levesque et al. 2010). 
The strange properties of the spectrum of GRB 090426 can 
therefore not simply be explained by a sightline effect. 

Considering the fact that this GRB exploded in what 
seems to be a galaxy merger, another scenario might explain 
the strange spectrum we observe. The lack of absorption 
from the normal ISM of the host suggests that the progen- 
itor did not reside in a star-forming region inside its host 
galaxy. A possible scenario is therefore that the progenitor 
was located in an isolated star-forming region outside the 
host itself, e.g. in something like a tidal tail created by the 
merging of the two galaxies. Thus, the sightline to the GRB 
would only intersect the material of the star-forming region, 
ionised by the GRB fiux and eventually piled up by the wind 
from massive Wolf-Rayet stars, and not the dense ISM of the 
galaxy. We note, however, that the afterglow position is, in 
projection, close to the center of the galaxy. 

Another strange feature about this burst was its rela- 
tively short duration paired with a soft high energy emission 
spectrum. Naturally, the question arises if the unusual high- 
energy properties and those of the afterglow spectrum could 
be connected. From Fig. [l] we see that all other low column 
density GRBs occupy the normal parameter space for long 
bursts concerning HR and duration, so there seems to be 
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no obvious correlation between HR or duration and low col- 
umn density. Since the material observed in the afterglow 
spectrum is relatively far away from the burst compared to 
the prompt emission, which is expected to be produced by 
internal shocks within the jet of the GRB, a connection is 
not to be expected. 

Concerning the progenitor of the burst, Antonelli et al. 
(2009) and Levesque et al. (2010) had already concluded 
that, despite its short duration, the burst was likely due 
to the collapse of a massive star. Arguments include the 
star-forming host galaxy, the high afterglow luminosity, con- 
sistency with the Amati relation, and the strength of the 
absorption lines. In past years, the classification of bursts 
purely on its duration or another single observational fact 
has also been questioned. Evidence for a coUapsar progen- 
itor has been observed even in the absence of a supernova 
signature (GRB 060505 and GRB 060614, see e.g. Fynbo et 
al. 2006b), which we cannot verify for GRB 090426 since the 
SN would be too faint to detect with current instruments, 
and a range of high redshift bursts have been detected with 
short intrinsic duration (see e.g. Lii et al. 2010). Suggestions 
for how to produce a coUapsar with such a short intrinsic 
duration include viewing the GRB from off-axis, such that 
we see only a small part of the jet (Lazzati et al. 2009), 
or the possibility that the central engine had been turned 
off before the jet reached the surface of the str, with the 
destruction wave catching up with the jet and leading to a 
shorter prompt duration (Mizuta & Aloy 2009). 

No other GRB observed so far has shown such extreme 
properties of the absorption lines, despite a number of spec- 
tra with relatively low column densities. However, none of 
them had time resolved spectral observations that would 
allow a line variability study (with the exception of GRB 
080310, a potentially similar case to GRB 090426 concern- 
ing the properties of the afterglow spectrum; see A. De Gia 
et al. in prep.). A larger number of rapid response observa- 
tions and time series of spectra using a high sensitivity and 
medium- to high-resolution spectrographs (in order to fit re- 
liable column densities) might allow us to study a few more 
of those rare events where we directly see the interaction of 
the GRB with its environment. Extremely rapid response, 
on the timescale of a few minutes, might even allow us to 
observe the material ejected by the progenitor itself. 
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